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In a class of non-singular cosmologies derived from higher-order corrections to the low-energy 
bosonic string action, we derive evolution equations for the most general cosmological scalar, vec- 
tor and tensor perturbations. In the large scale limit, the evolutions of both scalar and tensor 
perturbations are characterised by conserved quantities, the usual curvature perturbation in the 
uniform-field gauge and the tensor-type perturbed metric. The vector perturbation is not affected, 
being described by the conservation of the angular momentum of the fluid component in the absence 
of any additional dissipative process. For the scalar- and tensor-type perturbations, we show how, 
given a background evolution during kinetic driven inflation of the dilaton field, we can obtain the 
final power spectra generated from the vacuum quantum fluctuations of the metric and the dilaton 
field during the inflation. 
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I. INTRODUCTION 

Recent observations jjj of anisotropics in the cosmic 
microwave background radiation strongly support the 
paradigm of inflation as being the source of the primor- 
dial density fluctuations, producing an almost flat power 
spectrum [^j. The usual inflation models involve an ac- 
celerated period being driven by the energy density asso- 
ciated with the inflaton potential. However, there exist 
a class of models which lead to inflation simply because 
they are driven by the kinetic energy of a scalar field. 
There are a large class of such pole-like inflation mod- 
els: those derived from induced gravity scalar-tensor 
gravity or the pre-Big Bang scenario (PBB) j|] of 
string theory. Unfortunately, most of these models face 
the problem that they tend to yield blue spectra for both 
the scalar- and tensor-type perturbations with spectral 
indices ns — 4 and ut — 3, respectively, whereas the 
observationally supported Zel'dovich spectra correspond 
to ns — 1 and Ut — [H-H]. A possible resolution of 



this problem for the scalar perturbations was proposed 
in [lCfl , where fluctuations of the axion field present in 
the low-energy string action can generate the observed 
spectral index. However, it remains the case that the 
tensor spectrum is difficult to reconcile with this class of 
kinetic driven inflation models. 

Perhaps the biggest issue facing gravity models is how 
to tackle the initial curvature singularity. In the con- 
text of PBB, an interesting suggestion was made that 
the graceful exit problem could be resolved by including 
the quantum backreaction effect pTj |. A different ap- 
proach uses an expansion in terms of the inverse string 
tension (a') and coupling corrections. This has already 
met with some success [ ^2|Jl3| . Given that there exist a 
class of non-singular cosmologies based on these higher- 
order corrections, it is natural to investigate the effect 



of these correction terms on the evolution of primordial 
fluctuations. In fact the impact of potential higher-order 
curvature corrections on the evolution of the perturba- 
tions has recently been studied in various inflationary 
models derived from string theory [ ^4|JT5| . In particu- 
lar, the classical evolution in the presence of a Gauss- 
Bonnet coupling term was considered in JT^JT^]. Effects 
of a string theory motivated axion coupling term g(<j))RR 
in the Lagrangian, where RR = T) tlvaT R^J-^RaTtK, was 
investigated in jl^]. The general feature that appears 
to emerge is that the higher-order curvature corrections 
generally flatten the spectral distributions of primordial 
vacuum fluctuations. Hence, one may imagine a scenario 
where the observationally relevant perturbations for the 
large scale structures leave the Hubble radius during a 
highly-curved regime, leaving the power-spectrum nearly 
scale- invariant . 

The aim of this paper is to calculate the vacuum fluc- 
tuations arising out of the most general non-singular so- 
lutions formed to date from the first order curvature cor- 
rections, corrections which could arise in the context of 
the massless bosonic sector of the low-energy effective 
action of string theory. 

The paper is organised as follows. In §|n] we introduce 
a general action including possible curvature corrections 
up to fourth order in derivatives and we derive the cor- 
responding field equations. §111 is devoted to studying 



the classical evolution of three types of perturbations. 
We derive closed form equations for both scalar- and 
tensor-type perturbations and determine the correspond- 
ing large-scale exact solutions. We also show that the 
vector-type perturbations are described by the conser- 
vation of the angular momentum of the additional fluid 
in the absence of dissipative processes. In §1V, we dis- 



cuss the quantum generation of the scalar and tensor per- 
turbations, and derive the corresponding power-spectra. 
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The general results are applied in §|v| to the particular 
case of the pre-Big Bang scenario of string cosmology. 
Finally, we discuss our main results and future applica- 
tions in §fv|. 



II. THE EVOLUTION EQUATIONS 

To keep things as general as possible, we consider the 
following D-dimensional actiory 



S 



-V{4>) + i (c) + L T 



(1) 



where f(cj),R) is an algebraic function of a dimension- 
less scalar field <fi an d the scalar curvature R. u>(4>) and 
V(4>) are general algebraic functions of </>. Through the 
Lagrangian L^ c \ we allow for the inclusion of terms with 
even higher numbers of derivatives such as contracted 
quadratic products of the curvature tensor, whereas L m 
is the Lagrangian of additional matter fields (e.g. flu- 
ids, kinetic components, axions, moduli, etc.) with 
its associated energy-momentum tensor T^„ defined as 
\\/—<j[r ia/ &9nv = SiV^gLm)- By choosing the appropri- 
ate parameters, the general action eq. ([!]) includes Brans- 
Dicke theory, non-minimally coupled scalar field, induced 
gravity, i? 2 -gravity, etc JTsj] . For instance, Einstein grav- 
ity with a minimally coupled scalar field corresponds to 
the case / = R, w = 1, and = 0. 

In a string theory context (see || for a recent review 
and references therein), the low-energy effective action is 
obtained with / = e~*R, u = -e~^>, V = and = 0. 
Higher-order corrections take the form of an infinite se- 
ries expansion with expansion parameter a 1 — X 2 , where 
X s is the fundamental string length scale. For the sake 
of simplicity, we shall restrict ourselves to the simplest 
extension of the lowest-order gravitational action ensur- 
ing that the equations of motion remain second-order in 
the fields. In such a case, the most general Lagrangian 
density at the next to leading order in the Regge slope 
reads pjj : 

L [c) = -^a'A£(0) [ Cl R 2 GB +c 2 G^^, v 

+c 3 n^^ + c 4 (^0 ;M ) 2 ], (2) 

where G^ v = R^ v — ^g^R is the Einstein tensor and 
R 2 GB ee R^ pa R llvpa - AR^R^ + R 2 is the well-know 



We adopt the convention (—,+,...,+), R vp 



R 



8ttG = 



f) r M + 



and set our units such that ft 



Gauss-Bonnet combination ensuring the ghost-free char- 
acter of the theory. In fixing the coefficients c,-'s, we re- 
quire that the full action agrees with the three-graviton 
scattering amplitude |20| , and thus impose the con- 
straints, £ = — e~^, C3 = — [c 2 + 2(ci + C4)]/2 and 
ci = —1, working in units a' = 1. A is an additional 
parameter allowing for different species of string theo- 
ries, A = —1/4, —1/8 for the Bosonic and Heterotic string 
respectively, and A = for superstrings. The inclusion 
of such higher-order corrections for the background evo- 
lution has recently been investigated in the context of 
the pre-Big Bang scenario of string cosmology, leading to 
a number of promising non-singular cosmologies which 
smoothly interpolate between the growing to the decreas- 
ing curvature regime |T^Jl^| . 

By varying eq. (^) with respect to the metric and the 
scalar field we can derive the gravitational field equation 
and the equation of motion for the scalar field <fi, respec- 
tively: 



FGt = 



-S» (FR — f + 2V) 



V V 1 V 1 V 5 



- t { ; ] 



0. 



(3) 
(4) 



where F = J 1 ^ ^ and represent the contribu- 
tions derived from the next to leading order corrections 
given by eq. (§), 



T< c >£ ee -a'xl - 4 Cl [(R^ UT + RvaS? - S^Rar) e 
+c 2 [^R aT - R\ VT )rV - K<t>^ a 

-RZ^d>:a 



£ [2D0 0» v - 



+ 2 PWtv + H;A'n + Wh*^] 



2 



K = -5^R 2 GB + ±R\ VT R aT - 2RR£ 
T { , c) = a'xl c^RIb 



(5) 



(6) 



+2£ 



-Ci 



Equations (§)-(@) form a complete set of covariant gener- 
alised Einstein equations which include higher-order cor- 
rections. Hence they extend the domain of validity of any 
solutions into the highly-curved regimes included in eq. 
(||) . This system of equations provides the framework to 
study both the cosmological background and evolution 
of metric and field perturbations. Note that N£f corre- 
sponds to a quadratic expression which vanishes in four 
dimensions on account of the algebraic identities satisfied 
by the Riemann tensor pl| . 



III. PERTURBED FIELD EQUATIONS 

From now on, we concentrate on the four-dimensional 
case, and consider the metric of a spatially homogeneous 
and isotropic model with the most general perturbations: 



ds 2 = -a 2 (1 + 2a) dr] 2 - a 2 (/3 , + drjdx 1 
eff (1 + 2<p) + 2 1Aj + 2C m + 2Cij 



(8) 



dx l dx J , 



where a(t) is the cosmic scale factor with dt = adrj. Latin 
letters denote space indices. a(x,t), (3(x,t), <p(x,t) and 
7(x, t) characterise the scalar- type perturbation, £?j(x, i) 
and Cj(x, t) are transverse (B l ,. = = C 1 ^) and repre- 
sent the vector-type perturbation, whereas Cij(x, t) is 
transverse and tracefree (C^- = = C|), and corre- 
sponds to the tensor-type perturbation. Indices are based 

(3) 

on g\j as the metric, and a vertical bar indicates a co- 
variant derivative based on g}y. 

We decompose the energy-momentum tensor of the ad- 
ditional matter and the dilaton field into 



2?(x,t) 

<Kx,t) = 



4>(t) + S 



<5T^(x,<), 
<x,t). 



(9) 
(10) 



An overbar indicates a background order quantity and 
will be omitted unless necessary. The three types of per- 
turbations decouple from each other due to the symmetry 
in the background field equations and the fact that we are 
working to linear order in the perturbations. Thus, we 
can handle them individually and we find it convenient to 
separate their respective contributions in the perturbed 
energy-momentum tensor. We use the superscripts (s), 
(v) and (t) for the scalar, vector and tensor parts, respec- 



tively, such that STp = ST 



(s)fj. 



ST 



ST 



Using 



a normalised (u^u^ = —1) four-vector and its associated 
projection tensor = g^+u^u^, a covariant decompo- 
sition of the (imperfect fluid) energy-momentum tensor 
into fluid quantities is given by [p2[ 



T^v = pU^U u + ph^ + q^u v + q v u^ + 7r M „, 



(11) 



(7) where the background is assumed to be made up of a 
perfect fluid. Here, we have defined p = T^u^u" as 



the energy density and p = ^T^W as the pressure. 
The energy flux and anisotropic pressure are, respec- 
tively, = —T l , a u u h'J 1 and tt^ = To-i-h^hl, and satisfy 
q^u^ — = TTfiuU^ 1 , t^^v = TiVju. Convenient decomposi- 
tions for the fluid four-velocity and the energy flux are, 
respectively, it M = — a(l + a, Sut) and q^ = a(0,Sqi). The 
content of the energy-momentum tensor to linear order 
then becomes 



rpO 

1 o 



(p- 



tip), 

T? = -(p + p)S Ui 
T} = (p + 5p)S) + 



Sqi 



(12) 
(13) 
(14) 



In general, the decomposition of the flux in eq. ( |ll| ) sug- 
gests the apparition of additional degrees of freedom. For 
instance, we may choose to consider the normal frame 
Sui = or the energy frame Sqi = 0. In linear perturba- 
tion theory, however, Sui and Sqi always appear together 
in a frame independent combination of the from of eq. 
( |l3| ) . Hence, we can investigate the gauge transformation 
properties of such a combination in a frame independent 
way [ ^9| . For later convenience, we split eqs. ( p^| , pT[ ) into 
the three-types of perturbations as 



T 



■o _ 



(p + p) (v. 



T) EE ( P + S P ) 5) + n i3) ] + 7T w ; + 7T m } . 



(v)i 



Mi 



(15) 
(16) 



The background equations in § III A are presented con- 



sidering the general spatial curvature K of the back- 
groun d, whereas for the perturbation equations in §[IIB- 
III D we assume a fiat background (i.e. K — 0). For the 



scalar-type perturbation we set T ^ = 0, whereas for 
the vector and tensor type perturbations we keep the 
contribution from the fluid energy-momentum tensor. 
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A. Background field equations 



B. Scalar- type perturbations 



From eqs. (|3|)-(^) we obtain a set of equations de- 
scribing the background evolution of the homogeneous- 
isotropic FLRW world model in our generalised gravity 
model: 



H 2 



H 



K 



6F 



uj<\> 2 + RF - f + 2V - 6HF 



-2T ° - 2T (c) ° ,, 



K _ 1 / 
a 2 ~ 2F\ 



u>4> 2 + HF — F + T ° - -T- 



T (c)0 l T (c)i 

3 i/: 



(17) 



(18) 



4> + + ±- (u^ 2 - + 2V t4> + T ( ; ] ) = 0, (19) 

where we denote the Hubble parameter as H = a/a, the 
Ricci scalar R = 6(2H 2 + H + K / a 2 ) and the higher-order 
contributions are given by 



T 



[c)0 



-a'X 



-c 2 ^ 2 3H 2 



Vlr^H ( H 2 



T 



(c)t _ 



-3a'A<^ 4ci 



K 



3 r ^ 

K 



(20) 



^[2H + 2,H 2 



+ 2£H i^H + H 

A((j)H + 2^H 



K 
7a 



+3c 2 



+C30 



-2^,11 ( 2H + 3H 2 + J 



0£ + 3£(/> -6£(4>H+ 2<t)H + 3<pH' 



(22) 



The solutions of the system of equations (|17|)-(P2|) pro- 
vide the cosmological "graviton-scalar field" background 
in which we will study the propagation of scalar field and 
metric perturbations. Equations (p7|)~(^9[) are not inde- 
pendent. We can derive eq. ( p~8| ) from eqs. (|l7 19), by 
using T ° + H(3T§ - T?) = which follows from energy- 
momentum conservation of the additional fluids or fields. 



To investigate the scalar-type perturbations there are 
many different temporal gauge conditions available for us 
to use. Due to the homogeneity of the background three- 
space the spatial gauge transformation is trivial al- 
though [3 and 7 change under the spatial gauge transfor- 
mation, they always appear together in a spatially gauge- 
invariant combination x = a(P + a i), whereas a and (p 
are spatially gauge- invariant. Except for the synchronous 
gauge condition, which happens to be the most widely 
used one, all the other fundamental temporal gauge con- 
ditions fix the gauge degrees of freedom completely, thus 
a variable in such a gauge condition uniquely corresponds 
to a gauge-invariant combination of the variable con- 
cerned and the variable used in fixing the gauge con- 
dition, see eq. (|23|) and below. The most suitable gauge 
condition depends on the type of the problem and in gen- 
eral is not known a priori, hence the gauge-ready method 
proposed in ]23| , |T9] ] appears particularly convenient since 
it allows for a flexible use of the various fundamental 
gauge conditions. 

In handling the perturbations involved with the scalar 
field in Einstein or generalised gravity theories [without 
the term in eq. (0)], it is known that the field fluc- 
tuation in the uniform-curvature gauge or equivalently 
the curvature fluctuation in the uniform-field gauge (se e 
below) allow the simplest analysis [pif . Following p7| , 
we derive an equation for a gauge-invariant combination 
ipsd> defined as 



<P5<j> = <P ft 



H 



(23) 



In the uniform-field gauge, which takes 5<p = 0, the 
gauge-invariant quantity (ps^ is identified with <p. Sim- 
(21) ilarly 5(j) v is the same as 6(j> in the uniform-curvature 
gauge which sets (p = 0. The gauge invariant combina- 
tion (fis^ was first introduced by Lukash in ^5j,^6|. If 
we have a solution in one gauge condition, the other so- 
lutions in the same gauge as well as the ones in other 
gauges can be derived from it as linear combinations. In 
the present case, choosing the uniform-field gauge also 
implies <5£ — 0. We consider a case with F = F(<fi), thus 
we have SF = 0. Ignoring the additional matter, that is 
we do not include any other components other than the 



scalar (dilaton) field, i.e. T „ = 0, and following the 
same steps described above eq. (3) in [p7| , we can derive 
a closed form equation describing the classical evolution 
of the scalar-metric perturbation. Explicitly, we find 

^(^ (- ^)'-' W ^* = 0, (24) 



where 



Q 



,(£+o^r 

' 2F+Q ( h 3) 



Q 



(*) 



H 



(») 



4 



,w = 



= 1 



2F-\ 



-Qi 



Q 



LU(f) 2 



Q 



(») 



(25) 



The contribution arising from the higher-order correc- 
tions is summarised in the quantities 



>l s) = a'X 



-4 Ci £h 2 + 2c 2 ^ 2 H + c 3 C0 3 



^6 



= a'X 



-8a£H + c 2 & 2 



a'Xcj) 2 



-3c 2 £fP + 2c 3 



(i - m 



Q ( d s) ee a'Xtf [ - 2c 2 £ff - 2c 3 + (<j> - &H 

Qi s) = a'X [ - 16d£ff + 2c 2 ^> + 2£4> - 2&H 
-4c 3 ^ 3 



a'X 



(26) 



In units a' = 1, our result reproduces those previously 
obtained for the particular case c% = — 1, c 2 = Cs = C4 = 
and A = -1/4 Q. 



C. Vector-type perturbations 

Rotational perturbations do not enter the (0, 0) com- 
ponents of the field equations, but they do contribute to 
the (0, i) and components. In fact, our scalar field 
perturbation does not couple directly with these vector- 
type perturbations. To investigate the influence of the 
scalar field, we consider an additional fluid present in the 
system. We can introduce the vector-type fluid quantities 
of the additional matter as 

ST {V) ° =(p + p)v^Y} v \ 8T {v) ] ee tt^Y (v) ;, (27) 

where 5^ and are vector-type harmonic func- 

tions ^8|. is the velocity variable related to the 

vorticity and 71» is the anisotropic stress. In the field 
equations, the vector perturbations always appear in a 
gauge-invariant combination Bi + aCi = ^f( v ^Y^ . From 
the (0, i) and components of the gravitational field 
equation we can derive^]: 



Il_Q{v)y{v) = (p + p ) w (») i 



(28) 
(29) 



where we have introduced 



ee F - ia'A (8diH - c 2 & 



(30) 



If we ignore the anisotropic stress of the fluid, tt^ v \ eq. 
(^) implies the conservation of the angular-momentum 
of the fluid, as 



(p + p)-a-»W(x,*) =£(x). 



(31) 



Notice that this result is independent of the generalised 
nature of our gravity model. The evolution of vorticity 
depends on our generalised gravity indirectly through the 
background evolution. The generalised nature, however, 
appears directly in connecting the vorticity to the metric 
perturbation in eq. (^8|). Equation (p9|), in fact, follows 
from the 0-component of the conservation of the fluid 
part of the energy-momentum tensor {T®' ,u = 0), hence 
is naturally independent of the field equation. Conserva- 
tion of angular momentum forces the rotational pertur- 
bation to decay in an expanding medium, which renders 
it generally cosmologically uninteresting in the context 
of linear perturbations in expanding universe. 



D. Tensor-type perturbations 

From the component of the field equation we can 
derive the linearised equation for the gravitational wave 
perturbation 



a 

where 



<rrp(t)i 



(32) 



Q 



f -r.'X 



F-^a'x(8c^+c 2 ^ 2 



(33) 



8T ^ includes contributions to the tensor-type energy- 
momentum tensor, i.e., transverse-tracefree anisotropic 
stresses which can possibly arise if we include additional 
imperfect fluids. Since the gauge transformation does 
not affect the transverse-tracefree parts, Cj and 5T ^ 
are naturally gauge invariant. In the context of string 
cosmology, F = —ui = — £ = and V = 0, we note 
that one easily recovers the a 1 corrections obtained in 



IV. QUANTUM GENERATION AND 
EVOLUTION OF PERTURBATIONS 

A. Classical evolution 



t Compared with |28| we have = v c , 7r ( -"- ) ee pirijP and 



Assuming that there is no additional matter, T{? = 0, 
we have been able to derive a closed form equation for 
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the linearised classical evolution of both scalar and tensor 
metric perturbations. In a unified formalism [^7f , the 
dynamics of a perturbed variable $ is typically governed 
by the wave equation, 



1 (a»Q4)' -4* = Q, 



a 3 Q 



(34) 



which may be readily derived from the second-order per- 
turbed action 



s' 2 >s 



(35) 



where (Q, s, $) € {(Q( s \ ««, tp H ) ; (QW, s «, Cj)}. In- 
troducing the quantities ^ = with 2 = a^/Q, the 
linearised wave equation can be written for each 
Fourier mode, ipk = z$>k, in terms of the eigenstates of 
the Laplace-Beltrami operator, Aipk = —k 2 ipk- Explic- 
itly, eq. (34) becomes 



V4' 



sk 2 - V{rf) i\> k = 0, V(j]) 



(36) 



where ' = d/drj. The "pump" field z accounts for 
the parametric amplification of the metric fluctuations, 
whereas s may be interpreted as an effective shift in the 
frequency. We note that the wave equation ([36|) can be 
deduced from the second-order perturbed action 



ii' + — t 2 

z 



sipAip 



d 3 xdrj, 



(37) 



where the kinetic term is diagonal, emphasising the 
canonical nature of the variable tp. To linear order, the 
decomposition in Fourier modes implies that each comov- 
ing wavenumber k evolves independently of other comov- 
ing modes, satisfying eq. ( |36| ) for all times. In a Gen- 
eral Relativity context and in spatially flat FLRW man- 
ifold, we recall that eq. (^) reduces to the equation of 
a minimally coupled massless scalar field for tensor per- 
turbations |30f , whereas a time-dependent effective mass 
is present in the scalar perturbation case |31|. In the 
present context of generalised gravity theory, however, 
the wave equation differs from the Klein-Gordon equation 
for two reasons. On the one hand, both scalar and ten- 
sor perturbations are coupled not only to the background 
scale factor, but potentially also to the scalar field back- 
ground through the algebraic function F = df(4>, R) / dR. 
Hence, the growth of the comoving amplitude of metric 
perturbations rises because of the joint contribution of 
the metric and the scalar field background to the pump 
field z = z(a,<f>). On the other hand, the higher-order 
corrections act as an additional source for the paramet- 
ric amplification of the metric perturbations. These de- 
viations from the General Relativity case manifest them- 
selves through the unique function Q, which encompasses 
all sources of modification. 

In the case where the background evolution undergoes 
power-law-type inflation, the external potential responsi- 
ble for the parametric amplification process f32] behaves 



as V{rj) ~ \q\ 2 , where we have set z ~ I77I 7 . For each 



comoving wavenumber, the wave equation ( p6|) then re- 
duces to 



V4' + k 2 



±s 



v 2 - 1/4 



|fcry|2 



ipk = 0, 



(38) 



where we have introduced s = \s\ and v = i |1 — 2^]. At 
lowest-order, i.e. neglecting the frequency shift and cor- 
rections in the pump field, the difference between scalar 
and tensor metric perturbations is encoded in the ar- 
gument v G \v s \ v t }, reflecting the background dynamics 
through the exponent 7 G {7 S ; j t } of the pump field. The 
dynamical behaviour of the canonical variable, given by 
eq. (|38|), can be divided into two asymptotic regimes. On 
the one hand, the power-law behaviour of the external 
potential implies that it vanishes in the asymptotic past 
(\krj\ ^> 1). Hence, eq. ( |38| ) yields two oscillating solu- 
tions, corresponding to negative and positive frequency 
modes respectively. On the other hand, the monotonic 
growth of the effective potential of the perturbation en- 
sures that it will soon drive the evolution of the metric 
perturbation (parametric amplification process). In the 
large wavelength (\kr]\ <^ 1) limit, which corresponds to 
scales larger than the Hubble radius, eq. (38) yields the 
solution 



fcfcfox) = C fc (x) + ,D fc (x) / -f 
= C fe (x) + D fc (x)M 1 - 2 r, 



(39) 



where C k (x) and D k (x) are constants of integration. For 
7 < i , the metric perturbation approaches a constant 

for 77 — >• For 7 > | however, the second term with 

D k grows in time, with an additional logarithmic factor 
appearing at 7 = \. More generally, for both scalar 
and tensor metric perturbations in the large-scale limit, 
sk 2 -C z" /z, thus ignoring the Laplacian term in eq. (|3^), 
the exact solution is given by 



* fc (j7,x) =C fe (x)+D fe (x) 



dr)' 



(40) 



implying the conservation of $fc in the large-scale limit. 
Since the coefficient Cfe(x) does not depend explicitly 
on V (</>), f(4>,R) or w ((/)), the perturbations are con- 
served independently of general changes in the back- 
ground equation of state and remain a priori conserved 
even under changes of the underlying gravity. Since we 
haven't specified a cosmological scenario yet, these re- 
sults remain valid for any four-dimensional FLRW back- 
ground model which may be derived from the general 
action eq. ([!]). For modes with physical sizes of order of 
the Hubble radius, however, we stress that the conser- 
vation of the perturbed variable $fc may be delayed in 
a situation where the background evolution becomes in- 
creasingly sensitive to the effects of the higher-order cor- 
rections. Indeed, for both scalar and tensor metric per- 
turbations, the time-dependent variables Q and s may 



G 



eventually be altered by the presence of the curvature 
corrections. In that case, the source Q of the effective 
potential for the perturbation may no longer grow mono- 
tonically and may even decrease at the onset or during a 
high-curvature regime. In the context of standard infla- 
tion models, a similar effect has been associated in |33] l 
with an enhancement of the amplitude of scalar metric 
perturbations just after the scales cross the Hubble ra- 
dius in some special situations. But it remains from eq. 
( [40| ) that the evolution of the perturbed variables should 
be frozen out for scales far beyond the Hubble radius. 
Apparently, there is also no restriction on the sign of the 
frequency shift occurring when the curvature (a') correc- 
tions dominate the background dynamics: s may become 
negative or infinite depending on its particular form. One 
inter pre tation of the denominator of the frequency shift 
(eq. ( |25|) and eq. ( |33| ) for scalar and tensor metric pertur- 
bations, respectively) is that as the a' correction terms 
become comparable to the lowest-order terms, it marks 
the breakdown of the approximation we have adopted, in 
that we are entering a regime where we should include all 
the higher-order corrections. Although we expect these 
modifications to be dependent on the type of gravity cho- 
sen for the background evolution, we now proceed to in- 
vestigate further the evolution of the perturbation vari- 
ables inside the Hubble radius, for we are interested in 
obtaining the spectrum as the variables leave the Hubble 
radius during inflation. 



B. Quantum generation of vacuum fluctuations 

In the semiclassical approximation |34j] , the perturbed 
parts of the fields and metric are regarded as quantum 
mechanical operators, whilst the background parts are 
considered as classical. In the Heisenberg representation, 
where the quantum operators (denoted with an overhat) 
carry the time-dependence, we have for instance 

</>(x, t) = (j>(t) + S$(x, t), v?(x, i) -> (£(x, t), 



<P6<p = ¥> ' 



H 



= S$- jj<p, etc. (41) 



Tensor metric perturbations, leading to the formation 
of a cosmological stochastic background of gravitational 
waves, have two polarisation states, whereas scalar metric 
perturbations have a single component. For conciseness, 
we shall restrict the description of the quantum genera- 
tion to tensor-type perturbations; being without the po- 
larisation states the scalar-type perturbation is simpler, 
and can be read from the tensor-type analyses by ignoring 
the tensor and the polarisation indices. It proves conve- 
nient to introduce a decomposition based on the two po- 
larisation states which, for a flat three-space background, 
is given by @|g| 



$y(x, t) 



d 3 k A 



x, t; k) 



(42) 



d 3 k 
(2tt) 3 / 2 



ek(t)ai k elf (k) 



•hue. 



Here, h.c. denotes the Hermitian conjugate, I = +, x 
represent the two polarisation states and Q>eu(,t) is the 
mode function for the tensor metric perturbation. The 
(k) form a base of polarisation states which yields 



The annihilation and creation 



eg^e^^k) = 26m 
operators a^k and aj k of each polarisation state satisfy 
the standard commutation relations on constant time hy- 
persurfaces 



[6«c,o<'k'] = = [a] k ,aj, k ,] 
[a<k)<tyk'] = <$«'<5 3 (k - k'). 



(43) 



By projecting <f>y (x, t) on the basis of polarisation states, 
we easily obtain the mode expansion for each polarisation 
state of the gravitational wave: 



d A k 
(2tt) 3 / 2 
d 3 k 



$ ij (x,t;k)eW "'(k) 



(2tt) 3 / 2 



[e ik - x ^e k {t)a ek + h.c.] 



(44) 



The mode function $ek(t) is a complex solution of the 
classical mode evolution equation. Replacing the classi- 
cal variable <& with the Hilbert operator $ in the second- 
order perturbed action eq. ([35|) leads to an equation 
for the mode function <3^k which satisfies essentially the 
same form as the one obeyed by the classical variable, 
namely eq. ([34]). From the quadratic effective action in 
the perturbed variable, we then derive the momentum 
7r$ { canonically conjugate to <fv, 



7T$ f (x, t) = 



a 3 Q$ £ (x,i). 



(45) 



This implies that the equal-time commutation relation 
[$^(x,t),7r* / (x / ,t)] = i<5 3 (x-x') leads to 



$,(x,<),$,(x',i) 



x'). 



(46) 



Requiring agreement between eq. (|43|) and eq. (|46|), the 
Wronskian of the mode function ^^(t) must satisfy 



(47) 



where the index t means that the derivatives in the Wron- 
skian are with respect to the cosmic time coordinate. The 
power-spectrum based on the vacuum expectation value 
(a^ 1 0) ee 0, Vk) is then 



|$/(x + r,t)$«(x,t)|0)e 



- 4kr d 3 r 



^l**(*)| a . 



(48) 
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In other words, the power-spectrum corresponds to the 
Fourier transform of the two-point correlation function of 
the metric fluctuations. In a space without any preferred 
direction, the two polarisation states of the gravitational 
waves are expected to contribute equally. Hence, we have 



p<5„ ( k ' *) = 2 E v c t (M) = 2 J2 h d *® 



2tt 2 



(49) 



The power-spectrum of scalar-metric perturbations can 
be derived in a similar manner, although no polarisation 
decomposition is required p7|. We find 



k 3 



(50) 



We now turn our attention to the spectral distributions 
\C & (t)\ and \ips4,k(t)\. 

Depending on the sign of the frequency shift s — ±s 
(which we shall consider constant in a first approxima- 
tion) , eq. ( P8[) becomes either a Bessel equation (for posi- 
tive s) or a modified Bessel equation (for negative s). For 
positive s, the normalised mode function solution of eq. 
( |38| ) corresponds to a superposition of Hankel functions 
of the first and second kind, 



c 1 (k)H?\x)+c a (k)H<?>(xj\, (51) 



with x = vlfc|?7| from which the usual condition between 



the coefficients | C2 1 2 — | ci | 2 = 1 is obtained from eq. (47) 



In the large-scale limit, \/Ifc|?7| <C 1, we have, for v ^ 
and v = respectively, 



p l/2 

i 



,1/2 



i h s-» /2 r( 



/Q2tt aH\ri\ T(3/2) 

2^A5T fk_^ 3/2 
2tt 



2 



3/2-1/ 



Cl\/Q 



]n(V\s\k\v\). 



(52) 
(53) 



Here we have considered the conventional choice of vac- 
uum C2 = 1 and ci = 0, corresponding to positive fre- 
quency in the asymptotic flat space time for r\ — ► — oo. 
An additional y/2 factor appears for the gravitational 

1 /2 

wave power spectrum T 3 - which follows from proper 

Cij 

consideration of the two polarisation states; see eqs. 

©[§■ 

For the case of negative s, we have the modified Bessel 
functions l v (x) and K v {x) as the two independent solu- 
tions of eq. (pi), that is 



iVQ 



ci (k) I v {x) +c 2 {h)K v (x) 



(54) 



with x = V§k\r]\ and an unusual condition upon the coef- 
ficients c\c2 — c\c\ — i arising from eq. (|47j). In the large 
scale limit, the power spectrum exhibits similar depen- 
dence on the wave number as in eqs. (p2h-(|3|). However, 



modes evolving inside the Hubble radius face an instabil- 
ity since 



I v {x) 



K v {x) 



(55) 



which may lead to invalidate our assumption of small 
perturbations, hence the use of linear approach to cosmo- 
logical perturbations. This is the condition we alluded to 
earlier, as marking a breakdown in the underlying action, 
requiring us to replace it with even higher-order correc- 
tions to enable us to deal properly with a high curvature 
epoch. Fortunately, a number of features emerge in the 
weak coupling regime which allows us to proceed further. 



V. AN EXAMPLE: THE PRE-BIG BANG 
SCENARIO 

The obscrvationally relevant scales are thought to have 
exited the Hubble radius within about 60 e-folds before 
the end of inflation. We have in mind a scenario where 
inflation occurs in our generalised gravity. This is then 
followed by a smooth transition to the ordinary radiation 
dominated era of the standard FLRW model. In the large 
scale limit, the perturbation variable $ is conserved in- 
dependently of changes in the underlying gravity theory, 
and the power-spectra based on the quantum vacuum ex- 
pectation value can be identified at a later epoch as the 
classical power-spectra based on the classical volume av- 
erage. Thus, eqs. (|52])-(|53]) are now considered valid for 
the classical power-spectra and the spectral index of the 
scalar and tensor-type perturbations become |37fl 



ns — 1 = 3 — 2zA;, ut = 3 — 2v t . 



(56) 



Here, v s and v t are to be evaluated at the time of Hub- 
ble radius crossing during inflation, and will thus directly 
reflect the kinematics of the background when the per- 
turbations exit the Hubble radius. 

As an application of our general results, we may con- 
sider in particular the pre-Big Bang scenario of string 
cosmology M. In our formalism, it corresponds to the 
case / = e~vR, £ = = — uj and V = 0. Examples of 
regular backgrounds have recently been obtained by sup- 
plementing the low-energy effective action of string the- 
ory by the kind of higher-order corrections given in eq. 
(0) [|l2|,^3). Fig. H is a typical example of non-singular 
evolution: when the curvature scale is of order of the 
inverse string scale, the a' corrections may eventually 
stabilise the background into a de-Sitter like regime of 
constant Hubble parameter H ~ 0(\j 2 ) and linearly 
growing (in cosmic time) dilaton. Then, quantum loop 
corrections (based on the string coupling expansion) or 
a non-perturbative potential (yet to be determined) may 
trigger an exit to the FLRW radiation-dominated phase, 
with constant dilaton field. 

It has been known for some time that this cosmologi- 
cal model derived from string theory predicts a very blue 
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1 

0.9 
0.8 - 
0.7- 
0.6- 
0.5 
0.4 - 
0.3 - 
0.2- 



H 

1 /3 34>/S t 



-10 10 20 30 40 50 60 70 

FIG. 1. Here we reproduce (from |2S|| ) a non-singular 
evolution for the Hubble parameter H = a/a and for 0/3, 
as a function of the number of e- folds, N = In a. The 
low-energy, dilaton-driven phase takes place approximately 
for — oo < N < — 3. After a short transition, this initial period 
is followed by a string phase with nearly constant Hubble pa- 
rameter and linearly growing dilaton, for 2 < N < 55. After a 
successful exit triggered by loop corrections, the background 
evolution enters the FLRW radiation-dominated phase at 
N ~ 68. 



spectra for both the scalar- and tensor-type perturba- 
tions • Recalling the (conformal time) expression of 
the tree- level solutions §38], 



ofa) = H?)*' 1 "^, fa) = -V31n(-»,), 



(57) 



the exponents of the scalar and tensor pump fields coin 
cide, 7 S = 7t = 1/2, hence v s — v t — 0. Using eq. (56) 
we then obtain 



ns - 1 = n T = 3, 



(58) 



whereas the observationally favoured Harrison-Zel'dovich 
type spectra correspond to ng ~ 1 and tit — 0. Hence, 
the pole-like acceleration expansion (in the string frame) 
of the pre-Big Bang scenario cannot, in principle, be con- 
sidered as the source for the observed fluctuations in the 
CMBR and also the observed large scale structures. A 
possible resolution of this problem for the scalar pertur- 
bations was proposed in |1C|| , where fluctuations of the 
axion field present in the low-energy string action can 
generate the observed spectral index. There also remains 
a possibility that the high-curvature regime can last long 
enough so that all the relevant scales that we observe 
today actually exit the Hubble radius during this era. 

To determine whether a high-curvature regime may or 
may not flatten the spectral slope, we have to follow the 
evolution of the large frequency modes, associated with 
small physical wavelengths. We start by recalling that 
during the high-curvature regime, both the Hubble rate 
and the dilaton field become time-independent, H = H 
and = 0. Hence, for both scalar- and tensor-type per- 
turbations, the potential reduces to its tree-level (TL) 
form and we can express the pump field in the convenient 
form, z = const x ztl- Leaving aside the frequency shift, 



i.e. setting s = 1 in eq. (pq), it is already well known 
J39| that the spectral slope of the high-frequency modes, 
crossing the Hubble radius in the high curvature, stringy 
regime, and re-entering in the radiation era, is fully deter- 
mined by the fixed point values 0, H . Indeed, during such 
a string phase, the scale factor undergoes the usual de Sit- 
ter exponential expansion, while the logarithmic evolu- 
tion of the dilaton, in conformal time, is weighted by the 
ratio <f>/H, i.e. <j>{rj) (4>/H) log(-ry)+const Q. By 

introducing the convenient shifted variable = — 3H 

and the ratio $ = <fi/H, and referring the spectrum to 
a fixed point allowing a subsequent (loop catalysed) exit 
p3| , i.e. — 1 < i? < 0, one easily finds that the exponents 
of the pump fields satisfy < 7 S = -f t — \ (1 + < 1/2. 
Hence, the spectral indices for the high-frequency modes, 
leaving the Hubble radius during the de-Sitter like era 
and re-entering in the radiation dominance epoch is given 
by 



2 < ns — 1 = tit < 3. 



(59) 



Although still not satisfying the observational bounds, 
they provide an indication of a possible resolution for 
suitable higher-order corrections. 

In Fig. 0, we illustrate the evolution of the frequency 
shifts given by eq. ( p5~|) and eq. (|33| ) in the regular back- 
ground of Fig. ^. In the asymptotic past, the curva- 
ture corrections are negligible, i.e. a' — > and we have 
s ( s ) = §(*) = i At the onset of the high-curvature 
regime, the behaviour of the frequency shift for scalar 
and tensor-type perturbations is strongly affected by the 
a' corrections, is attracted to a positive value while 
decreases and becomes even negative. Then, when 
the universe reaches the FLRW radiation-dominated era 
with a constant dilaton field, the higher-order corrections 
identically vanish and we recover — = 1. 

For tensor-type perturbations, the influence of the (all- 
time positive) frequency shift has been recently investi- 
gated in p9]| . Its contribution has been found to amount 
to an overall rescaling, by a numerical factor of order 
unity, of the total energy density of the background, but 
does not alter the spectral index of the high-frequency 
modes in eq. (|5|). For scalar- type perturbation, however, 
the situation is more complex, due to the sign change of 
the effective frequency shift during the high-curvature 
regime. For a comoving mode evolving inside the Hub- 
ble radius at the onset of the high-curvature phase, the 
sign change is a priori responsible for turning the oscil- 
latory behaviour of the perturbation variable into expo- 
nential growth or exponential suppression. In the latter 
case, this may imply an abrupt end-point in the spectral 
distribution, where the highest frequency mode to get 
amplified corresponds to the one which left the Hubble 
radius just before the onset of the string high-curvature 
phase. In the former case, on the contrary, the growth 
of the amplitude of the perturbation may be exponential 
until it exits the Hubble radius. Neither case seems sat- 
isfactory, and the change of sign is really an indication 
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-10 10 20 30 40 50 60 70 

FIG. 2. Here we compare the frequency shifts for scalar 
and tensor-type perturbations, s' s ' and s'*' respectively, in 
the regular background of Fig. hi The frequency shift for ten- 
sor-type perturbations is found to be always positive, whereas 
a 1 corrections are responsible for the sign change of the fre- 
quency shift of scalar-type perturbations. The latter may 
imply an exponential instability for the perturbation am- 
plitude of the comoving mode whose physical wavelength is 
smaller than the size of the Hubble radius during the stringy 
high-curvature regime. 



that the scalar fluctuations should be considered in the 
light of the full string action in this regime. For example, 
including higher-order curvature contributions to probe 
the highly-curved regime, we have focused on the next 
to leading order in the infinite series of the inverse string 
tension expansion. Although this represents a significant 
improvement compared to previous studies, it still places 
a limit to our ability to analyse the results. It is natural 
to ask how sensitive our results are with respect to the 
inclusion of even higher-order derivative terms, such as 
R 3 or R^.pRV'f . To estimate the impact of such terms, 
while conserving equations of motion to second-order in 
the fields, we have looked at the relative contribution of 
(i9 M 0) 4 and (d^cj)) 6 . Indeed this very preliminary analysis 
indicates that the frequency shift may no longer be nega- 
tive by including the infinite series expansion. The tensor 
contributions are not affected by this kind of higher order 
terms. 



ing ground work for future applications to the quantum 
generation processes. In particular we have calculated 
the final spectra of the scalar- and tensor-type perturba- 
tions generated from the vacuum quantum fluctuations 
of the metric and the field, given the non-singular back- 
ground evolution during the kinetic driven inflation. 

Our results are valid for the general action of eq. ( 1 ) , 
and for the general perturbed FLRW metric in eq. (£). 
Moreover our results up to §IV describing the classical 
evolution of the background and the perturbations are 
generally valid without assuming any cosmological sce- 
nario, like the pre-Big Bang or other types of inflation 
models. In order to calculate the initial conditions for 
the seed structures generated from the quantum fluctu- 
ations, we need to specify the specific cosmological sce- 
nario for the background. As an example, we have ap- 
plied our general results to the case of the pre-Big Bang 
scenario of string cosmology. Our results confirm previ- 
ous expectations, that the low frequency modes, crossing 
the Hubble radius in the low-curvature regime, are unaf- 
fected by higher-order corrections. However, the higher- 
order corrections generally reduces the spectral indices 
for the high-frequency modes which are leaving the Hub- 
ble radius during the de-Sitter like era. Although still not 
satisfying the observational bounds, they provide an in- 
dication of a possible resolution for suitable higher-order 
corrections. 
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VI. CONCLUSION 

In this paper, we have for the first time presented a 
complete set of equations describing the classical evo- 
lution of scalar, vector and tensor perturbations in the 
context of a non-singular cosmology arising out of higher- 
order corrections to the low-energy string action. From 
these we have obtained the large-scale exact solutions 
for both the scalar and tensor perturbations, the asso- 
ciated large-scale conserved quantities and shown how 
there is general conservation of the angular momentum 
for the vector perturbations in the absence of dissipa- 
tive processes. The perturbed actions for the scalar- and 
tensor-type perturbations have been obtained, thus lay- 
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